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A review of the approximations in any time-depedendent density functional calculation of excitation 
energies is given. The single-pole approximation for the susceptibility is used to understand errors 
in popular approximations for the exchange-correlation kernel. A new hybrid of exact exchange and 
adiabatic local density approximation is proposed and tested on the He and Be atoms. 



I. INTRODUCTION 

Ground-state density functional theory is well- 
established as an inexpensive alternative to traditional 
ab initio quantum chemical methods Q. Now time- 
dependent density functional theory (TDDFT) is rapidly 
emerging as an inexpensive accurate method for the 
calculation of electronic excitation energies in quantum 
chemistry 1^,^. Calculation of dynamic response prop- 
erties using TDDFT has a long history, since the pio- 
neering work of Zangwill and Sovcn |^,|^. It is only 
relatively recently that attention has been focussed on 
the direct extraction of excitation energies [|6|-|ll|]. Al- 
ready, this method has been implemented in several 
quantum chemistry packages, such as dcMon Tur- 
bomole ADF and QCHEM Impor- 

tant calculations include the calculation of excited-state 
crossings in formaldehyde p^ , excitations with signif- 
icant doubly-excited character |l^], the photospectrum 
of chlorophyll A jl^, and the response of 2-D quantum 
strips 

How are excitation energies calculated using TDDFT? 
First, a self-consistent ground-state Kohn-Sham calcu- 
lation is performed, using some approximation for the 
exchange-correlation energy E^c, such as B3LYP or 
FEE [Q. This yields a set of Kohn-Sham eigenvalues 
ti and orbitals (pi. Even with the exact ground-state 
energy functional and potential, these eigenvalues are 
in general not the true excitations of the system, but 
are closely related. In a second step, the central equa- 
tion of TDDFT response theory is solved, which extracts 
the true linear response function from its Kohn-Sham 
counterpart. This equation includes a second unknown 
functional, the exchange-correlation kernel /xc(r,r';w), 
which is the Fourier transform of the functional deriva- 
tive of the time- dependent exchange-correlation poten- 
tial. The poles of the exact response function are shifted 
from those of the KS function, and occur at the true exci- 
tations of the system. These steps are typically repeated 
for several nuclear positions. 



The success of any density functional method, however, 
depends on the quality of the approximate functionals 
employed. The above calculation requires two distinct 
density functional approximations: one for the ground- 
state energy, which implies a corresponding approxi- 
mation for the exchange-correlation potential Uxc(r) = 
6Exc/Sp{r), and a second for the exchange-correlation 
kernel. Most calculations now appearing in the chemical 
literature use the adiabatic local density approximation 
(ALDA) for /xc- Adiabatic implies that the frequency- 
dependence of /xc is ignored, and its = value used, 
and LDA implies /xc(r, r') = 5wxc(r)/5p(r'), where v^c 
is the LDA potential from ground-state LDA calcula- 
tions. 

The more drastic of the two approximations is that for 
the ground-state. Very often, functionals which yield ac- 
curate ground-state energies have very poor-looking po- 
tentials. (How this can happen can be understood by 
considering the virial theorem, which relates energies to 
potentials These poor potentials then have in turn 

badly behaved virtual orbitals. In particular, local and 
semi- local approximations, i.e., generalized gradient ap- 
proximations, fail to capture the correct asymptotic be- 
havior of the potential, and many virtual states are not 
even bound. Even hybrid functionals, such as B3LYP 
(]l9| and PBEO [Q, do not much improve the asymp- 
totic behavior, since they only mix a fraction of exact ex- 
change. This restricts calculations within these approx- 
imations to only low-lying excitations. This difficulty 
is most pronounced in atoms, becomes smaller for bigger 
molecules, and is irrelevant for bulk solids. Recently [p3| , 
we have shown that even with approximations which are 
free from self-interaction error, such as exact exchange, 
or self-interaction corrected LDA, and which therefore 
reproduce the dominant part of the asymptotic decay of 
the potential, inaccuracies in the Kohn-Sham eigenvalues 
(mostly due to the incorrect position of the highest oc- 
cupied level) dominate over any errors introduced in the 
second step, namely the correction of KS levels to the 
true levels. 
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The present work studies the accuracy of approxima- 
tions to /xc alone. This is because, for low-lying states, 
errors in /xc must be disentangled from errors in the KS 
energy levels, and also because many people are work- 
ing to improve approximations to the exact ground-state 
potential p^-[2^, which we hope will ultimately reduce 
those errors discussed above. Hence all our calculations 
are performed using the exact Kohn-Sham potentials of 
the He and Be atoms, for which we thank Cyrus Umrigar 
p7t . We study a variety of approximations to /xc , all adi- 
abatic. We focus especially on the spin-decomposition of 
such approximations, which determines the relative po- 
sitions of singlet and triplet excitations in TDDFT. We 
use the single-pole approximation for the susceptibility to 
directly relate errors in /xc to errors in excitation ener- 
gies. We find that exact exchange works well for parallel 
spins, which determines the mean energy of the singlet 
and triplet, while the antiparallel contribution to /xc de- 
termines their splitting, which is well-approximated in 
ALDA. With this insight, we construct a hybrid of ex- 
act exchange and ALDA, which greatly improves results 
for He, and moderately improves them for Be. We use 
atomic units (e^ = Ti = rUe = 1) throughout, except in 
Fig. 1. 

II. METHODOLOGY 

The basic response equation of TDDFT has the same 
form as that of time-dependent Hartree-Fock theory, or 
of the Random Phase Approximation, i.e., 

X^^'(r,r';u;)=xr'(r,r';c^) 

+ E / J d'r'" xr"ir,r";Lo) 

<j" a'" 

/,C""'(r",r"';^)x^"'^'(r"',r';a;) (1) 

where x is the exact frequency-dependent susceptibility 
of the system, while Xs is its Kohn-Sham analog, and 
jli^{v,v') = l/|r-r'|-|-/^^'(r,r'). The cr indices denote 
spin, i.e., g =| or \. By various means, the poles of x as 
a function of lo can be found . Our method for finding 
these poles is to consider only the discrete poles, i.e., 
those corresponding to bound states. In the particular 
case of a frequency-independent model for /xc, we can 
show 1^ that these poles occur at the eigenvalues of the 
matrix 



qa 



+ aq.,> d?r /dV$* (r) x 



(2) 



/„-x"c(r,r')$9V'(r'). 

For notational brevity, we have used double indices q = 
(j, k) to characterize the excitation energy LOq^ = ija—^ka- 
of the single-particle transition {ka —> ja). Conse- 
quently, we set Q!go- := fka — fja, where fka- is the occupa- 
tion number of that orbital, and <i>qcr(r) = (pl.^{r)ipjcr{r). 



While we do include sufficient bound-state poles to con- 
verge to an accurate result, our method does neglect con- 
tinuum contributions, and this effect will be discussed in 
the next section. 

All approximations we study for /xc are adiabatic. The 
most ubiquitous is ALDA (or more precisley, the adia- 
batic local spin density approximation) in which 



/•crcr'ALDA 
J xc 



(r,r')=^(r-r') 



duaduo 



"T('"),"i(r') 



(3) 



Note that this leads to a completely short-ranged approx- 
imation to /xc- Similarly, any adiabatic GGA approxi- 
mation leads to an approximate /xc which is almost as 
short-ranged. 

A second distinct approximation to /xc is in terms of 
its (usually) dominant exchange contribution. A highly 
accurate approximation to the exact exchange-only equa- 
tions of ground-state density functional theory (the op- 
timized effective potential equations) was introduced by 
Krieger, Li, and lafrate |2^. This approximation has 
been extended to the time-dependent case H: 



n„{r)\r - r'\n„{r') 



(4) 



This is exact for one electron, and for (spin-unpolarized) 
two-electron exchange. Note that this approximation has 
a long-ranged contribution, which can cancel exactly the 
direct hartree contribution to the matrix M in Eq. (H). 

A third approximation which we tried is the self- 
interaction corrected (SIC) ALDA, which is simply the 
second functional derivative of the SIC-LDA energy: 



E: 



where n,- 



iSIC~LDA 



{(Pial"^ is the density of a single-orbital, and 
Ex[nia] is the exact Hartree self- interaction energy of 
that orbital. This approximation should improve over 
ALDA, in avoiding spurious self-interaction errors, and 
over just exchange, by including some correlation. 



III. DATA 

In this section, we report calculations for the He and Be 
atoms using the exact ground-state Kohn-Sham poten- 
tials, the three approximations to the kernel mentioned 
in the previous section, and including many bound-state 
poles in Eq. (^), but neglecting the continuum. The 
technical details are given in Ref. |^^. Table I lists the 
results, which are compared with a highly accurate non- 
relativistic variational calculations ||2^,^ In each sym- 
metry class (s, p, and d), up to 38 virtual states were 
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TABLE I. Singlet/triplet excitation energies in the helium 
and beryllium atoms, calculated from the exact Kohn-Sham 
potential by using approximate xc kernels (in millihartrees) , 
and using the lowest 34 unoccupied orbitals of s and p sym- 
metry for He, and the lowest 38 unoccupied orbitals of s, p, 
and d symmetry for Be. Exact values from Ref. for He 
and from Ref. g for Be. 



Singlet/triplet shifts 






ALDA 


X 


SIC 


hybrid 


exact 




Transitions from the Is state in 


He atom 




2s 


746.0 


22/- 11 


20/-25 


19/-16 


14/-19 


12/-18 


3s 


839.2 


6.9/-2.4 


5.8/-4.9 


5.6/-3.6 


5.1/-4.6 


3.3/-4.2 


4s 


868.8 


3.1/-0.9 


2.5/-1.7 


2.4/-1.3 


2.4/-2.0 


1.3/-1.6 


5s 


881.9 


1.6/-0.4 


1.3/-0.8 


1.3/-0.6 


1.3/-1.0 


0.6/-0.8 


6s 


888.8 


1.0/-0.3 


0.8/-0.5 


0.7/-0.4 


0.8/-0.7 


0.4/-0.5 


2p 


777.2 


-0.8/-7.4 


7.2/-8.4 


6.1/ 0.2 


2.7/-3.4 


2.7/-6.6 


3p 


847.6 


0.7/-1.9 


2.5/-2.3 


2.2/-0.5 


1.4/-1.3 


1.0/-2.0 


4p 


872.2 


0.4/-0.7 


1.1/-0.9 


1.0/-0.2 


0.6/-0.5 


0.5/-0.8 


5p 


883.6 


0.2/-0.4 


0.6/-0.5 


0.5/-0.1 


0.4/-0.3 


0.2/-0.4 


6p 


889.8 


0.1/-0.3 


0.3/-0.3 


0.3/-0.1 


0.2/-0.2 


0.1/-0.3 


err 


57 


32 


31 


29 
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Transitions from the 2s state in 


Be atom 




oS 


244.4 


7.1/-5.7 


10.9/-10.6 10.3/-1.4 


6.6/-4.6 


4.7/-7.1 


4s 


295.9 


2.5/-1.6 


3.6/-2.5 


3.5/-0.6 


2.6/-1.6 


1.4/-2.0 


5s 


315.3 


1.1/-0.7 


1.7/-1.0 


1.6/-0.3 


1.2/-0.7 


0.6/-0.9 


6s 


324.7 


0.6/-0.4 


0.9/-0.5 


0.9/-0.2 


0.7/-0.4 


0.3/-0.5 


2p 


132.7 


56/-42 


55/-133 


53/-53 


10/-88 


61/-32 


3p 


269.4 


2.0/-4.3 


6.4/-4.2 


5.6/ 1.1 


4.2/-2.8 


4.8/-1.5 


4p 


304.6 


0.3/-1.4 


2.1/-1.2 


1.9/ 0.3 


1.3/-0.8 


1.7/-4.1 


5p 


319.3 


0.1/-0.6 


1.0/-0.5 


0.9/ 0.1 


0.6/-0.2 


0.2/ 0.0 


6p 


326.9 


0.0/-0.4 


0.5/-0.3 


0.4/ 0.0 


0.3/-0.2 


0.1/-0.1 


3d 


283.3 


-5.4/-2.8 


1.8/-2.0 


0.9/ 3.2 


-1.4/ 1.2 10.3/-0.6 


4d 


309.8 


-1.4/-1.1 


0.8/-0.9 


0.5/ 0.9 


-0.2/ 0.6 


3.6/-0.2 


err 


138 


56 


144 


73 


136 




err' 


45 


41 


37 


44 


29 





calculated. The errors reported are absolute deviations 
from the exact values. The second error under the Be 
atom excludes the 2s — > 2p transition, for reasons dis- 
cussed in the next section. 

The effect of neglecting continuum states in these cal- 
culations has been investigated by van Gisbergen et al. 

who performed ALDA calculations from the exact 
Kohn-Sham potential in a localized basis set. These cal- 
culations were done including first only bound states, 
yielding results identical to those presented here, and 
then including all positive energy orbitals allowed by 
their basis set. They found significant improvement in 
He singlet- singlet excitations, especially for Is — > 2s and 
Is — !■ 3s. Other excitations barely changed. Assuming 
inclusion of the continuum affects results with other ap- 
proximate kernels similarly, these results do not change 
the basic reasoning and conclusions presented below, but 
suggest that calculations including the continuum may 
prove to be more accurate than those presented here. 



IV. SINGLE-POLE ANALYSIS 

The simplest truncation of the eigenvalue equation (j^) 
for the excitation energies is to ignore all coupling be- 
tween poles, except that between a singlet-triplet pair. 
This is equivalent to setting ((jI/hxcI?') to zero, for q ^ q' . 
(We have dropped the spin-index on these contributions, 
since we deal only with closed shell systems). Then the 
eigenvalue problem reduces to a simple 2x2 problem, 
with solutions 



r!+ =cj, + 2K(9|/Hxck) 
r!- =c^, + 25R(9|A/xc|<7), 



(6) 



where 



/hxc — ^ ^ /hxc — Ij. _ +2 ^'^^'^ ^ ^^'^ ' 

(7(7' 



(7) 



Thus /hxc is the spin-summed contribution, which con- 
tributes to Xj the spin-summed susceptibility, and there- 
fore gives rise to the singlet level, while A/xc is the spin- 
flip contribution, also called /i^Gxc in the theory of the 
frequency-dependent magnetization density ]3H |. Thus 
even within the SPA, the KS degeneracy between sin- 
glets and triplets is broken, and we identify with the 
triplet. In Table II, we report results within the single- 
pole approximation. 

At this point, we notice the very strong shift in the Be 
2s — > 2p transition. This is due to the small magnitude 
of its transition energy, so that the pole of the 2p energy 
is very close to the pole of the 2s energy. Thus the single 
pole approximation is not expected to work well for this 
case, and it should be excluded from general statements 
based on the SPA. 



A. Why are Kohn-Sham excitation energies so good? 

We see throughout the data that the Kohn-Sham 
eigenvalues are always inbetween the exact singlet and 
triplet energy levels. The splitting is much larger for Be 
than for He, but this obervation is true in both cases. 
It has already been made by Filippi et al |^|, and ex- 
plained in terms of quasi-particle amplitudes |3^. Here, 
we use the single-pole approximation to analyze this re- 
sult in terms of the known behavior of density function- 
als. From Eq. (6) we find that the mean energy is given 
by 



1 



/IJ(r,r'|g), 



(8) 
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TABLE II. Same as Table I, but within the single pole 
approximation. 









Singlct/triplct shifts 






exact 


ALDA 


X 


SIC 


hybrid 




Transitions from the Is state 


in He atom 




2s 


746.0 


25.8/-10.3 


22.7/-22.8 : 


21.4/-14.7 


18.0/-18.1 


3s 


839.2 


6.6/-2.6 


5.6/-5.5 


5.3/-3.9 


4.7/-4.5 


4s 


868.8 


2.6/-1.0 


2.2/-2.1 


2.1/-1.5 


1.8/-1.8 


5s 


881.9 


1.3/-0.5 


1.1/-1.1 


1.0/-0.8 


0.9/-0.9 


6s 


888.8 


0.7/-0.3 


0.6/-0.6 


0.6/-0.4 


0.5/-0.5 


2p 


777.2 


-0.8/-7.0 


7.8/-7.9 


6.4/ 0.2 


3.1/-3.1 


3p 


847.6 


0.7/-2.0 


2.4/-2.3 


2.1/-0.5 


1.4/-1.3 


4p 


872.2 


0.4/-0.8 


1.0/-1.0 


0.9/-0.2 


0.6/-0.6 


5p 


883.6 


0.2/-0.4 


0.5/-0.5 


0.5/-0.2 


0.3/-0.3 


6p 


889.8 


0.1/-0.3 


0.3/-0.3 


0.2/-0.1 


0.2/-0.2 


err 


57 


34 


32 


31 
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Transitions from the 2s state 


in Be atom 




3s 


244.4 


8.2/-5.4 


13.0/-9.5 


12.2/-1.3 


8.6/-5.1 


4s 


295.9 


2.4/-1.6 


3.5/-2.7 


3.4/-0.6 


2.4/-1.6 


5s 


315.3 


1.0/-0.7 


1.5/-1.1 


1.4/-0.3 


1.1/-0.7 


6s 


324.7 


0.5/-0.4 


0.7/-0.6 


0.7/-0.2 


0.5/-0.4 


2p 


132.7 


75/-34 


71/-65 


67/-41 


58/-52 


3p 


269.4 


-0.4/-4.3 


5.4/-4.7 


4.0/ U.D 


2.3/-1.6 


4p 


304.6 


-0.1/-1.5 


1.7/-1.5 


1.5/ 0.2 


0.8/-0.6 


5p 


319.3 


-0.1/-0.7 


0.8/-0.6 


0.7/ 0.1 


0.4/-0.2 


6p 


326.9 


-0.1/-0.4 


0.4/-0.4 


0.3/ 0.0 


0.1/-0.2 


3d 


283.3 


-5.0/-2.6 


1.8/-1.9 


1.0/ 3.4 


-1.2/ 1.1 


M 


309.8 


-1.4/-1.1 


0.8/ 0.0 


0.6/ 0.8 


0.3/ 0.6 


err 


138 


61 


79 


60 


55 


err' 


45 


45 


36 


44 


33 



5?(gl/c"k) 



(2el) 



(12) 



for two electrons in the single-pole approximation. It is 
well-known (see, e.g., Ref. |Q) that for ground-state en- 
ergies, parallel correlation is much weaker than antiparal- 
lel, since antiparallel electrons are not kept apart by the 
exchange interaction. Thus this ratio is expected to stay 
well less than 1, as it does for all our He excitations. The 
effect of the single-pole approximation on this conclusion 
can be judged by studying the shift in the mean for the 
X results in Table I. 

For Be, and any system with more than two electrons, 
there is still a good deal of cancellation of the exchange 
contribution with the direct contribution, but this can- 
cellation is no longer exact. This can be seen in the Be 
results for X in Table II. By studying the form of /x 
given in Eq. (4), we expect this remnant exchange con- 
tribution to be of order 0{{N — 2)/N^) for unpolarized 
systems. Thus in exchange-dominated (i.e., high density 
or weakly correlated) systems, the direct Coulomb term 
in the denominator will be larger than any remnant ex- 
change term in the numerator. On the other hand, in 
low-density or strongly correlated systems, if antiparallel 
correlation continues to dominate over parallel correla- 
tion, this ratio will still be less than one. We conclude 
that the Kohn-Sham levels will usually be close to the 
true excitations (of single-particle nature). 



while the energy splitting is given by 

An, = in+ ~n-)^2mj^^^+flHQ), (9) 

since there is no exchange contribution to antiparallel 
/xc- If we further define Sflq — Qq — ujq as the deviation 
of the mean energy from the Kohn-Sham level, we see 
that 

\Snq\<Anq/2 (10) 

must be satisfied for the Kohn-Sham level to lie in be- 
tween the singlet and triplet levels. Within the single- 
pole approximation, we have a very simple expression for 
the ratio of these two: 

2dnq _ ^{qlflUl) 

Consider first the He atom. For two electrons, /l^ = 
— l/|r — r'l, exactly cancelling the Hartree term in Eq. 
(^, leaving only the parallel correlation contribution. 
(This is reflected in the X column on Table II, where 
the upshift of the singlet is equal to the downshift of the 
triplet.) Thus we find 



B. Relation of exact exchange to Gorling-Levy 
perturbation theory 

Both time-dependent DFT and Gorling-Levy |^ per- 
turbation theory are formally exact methods for extract- 
ing electronic excitation energies in density functional 
theory. In this section, we consider the expansion of the 
excitation energies in powers of the adiabatic coupling 
constant A to first order. This procedure should give 
identical results in both theories. Recently, Filippi et al 
p2| have performed GL first order calculations for the 
He atom, using the exact Kohn-Sham potential. Their 
results are numerically identical to ours, using the exact 
exchange kernel, but only within the single-pole approx- 
imation, as given by Table II. Results calculated with 
the full (i.e., many poles) scheme, i.e., in Table I, differ 
slightly from theirs. This produces a paradox, in which 
the easier SPA is more accurate (apparently exact), while 
the more sophisticated treatment introduces errors. 

The resolution of this paradox can be seen most easily 
in Eq. (1), the RPA-type equation for the susceptibility. 
Insertion of /x alone (linear in A) into these equations 
will lead to all powers of A being present in the solution, 
since it is a self-consistent integral equation. This is most 
easily seen by iterating the equations. A simple way to 
recover the exact first-order GL result is by solving the 
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FIG. 1. Excitation energies from the ground-state of He, 
including the orbital energies of the exact Kohn-Sham po- 
tential, the time-dependent spin-split correction calculated 
within four different functional approximations, and exper- 
imental numbers ISSl. 



equations with A/x, and making A very small, to find the 
linear contribution to the change in excitation energies. 
Insertion of A = 1 into this result will yield the exact GL 
result. 

A far simpler method is to use the single-pole approx- 
imation. To see why this works, consider Eq. (2), our 
matrix whose eigenvalues are at the excitation energies. 
Since /hxc is at least first-order in A, all off-diagonal con- 
tributions are of 0(A). Thus to lowest order, the di- 
agonal dominates, and the off-diagonal corrections con- 
tribute O(A^) corrections. Retaining only diagonal con- 
tributions, i.e., the single-pole approximation, yields the 
exact result to first order in A. A detailed functional 
derivation of this result has recently been given by Gonze 
and Scheffler [57|. 



C. A new hybrid functional for the kernel 

To illustrate the importance of understanding the ori- 
gin of errors in density functionals, we use the insight 
gained within the SPA in the previous sections to con- 
struct a hybrid functional for f^c- We then apply this 
both within SPA and in the full calculation. Our original 
idea, as mentioned in section II, was to use a simple self- 
interaction correction to produce a better approximation 
than either exact exchange or ALDA, but the SIC results 



of Tables I and II show this does not happen. 

Consider Fig. 1, which illustrates the positions of en- 
ergy levels in the different schemes for the He atom. Our 
first step is to consider the mean energies. As pointed out 
in subsection A, these are determined by the parallel con- 
tributions to /xc and the Hartree term. We have already 
seen how these two terms cancel exactly at the exchange 
level, so that the mean energy is very good in such a cal- 
culation. In ALDA, the cancellation (exact for He) of the 
exchange contributions is lost. This can be seen in the 
large shifts in the ALDA mean energies in Table II and in 
Fig. 1. Furthermore, the remaining small parallel-spin 
correlation contribution can be expected to be grossly 
overestimated by ALDA, since LDA ground-state corre- 
lation energies are usually too large by a factor of 2 or 
3. Unfortunately, even SIC- ALDA is not exact for two- 
electron exchange, and it also suffers from a poor mean 
energy. Thus we recommend using only exact exchange 
for the parallel-spin contribution. 

On the other hand, the splitting is determined solely 
by anti-parallel correlation contributions, and the direct 
term. Thus, an exact exchange treatment misses en- 
tirely the significant anti-parallel correlation contribu- 
tion. This error is highlighted by the far too large split- 
tings in the exchange results in Tables I and II, implying 
significant cancellation between the direct and antiparal- 
lel correlation contributions. So here we advocate use of 
ALDA. Since the splitting depends on anti-parallel con- 
tributions to /xc, but the SIC correction only applies 
to one spin at a time, SIC- ALDA has exactly the same 
splittings. Our recommended hybrid is therefore 



fTT = fTT 



fTi ^ fTi ALDA 
J xc J xc 



(13) 
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Results shown in Tables I and II indicate that this hybrid 
decreases almost all errors over either exact exchange or 
ALDA. On average, the decrease is by about a factor of 
3 for the He atom, but much less for Be (about 40%). 
Notable exceptions are the triplet transitions to p states 
in He and in Be, where the error is increased. 

V. CONCLUSIONS 

We have shown how the results of TDDFT with ap- 
proximate exchange-correlation kernel functionals may 
be understood in terms the well-known behavior of the 
ground-state functionals from which they are derived. 
We have shown in a simple case how a more accurate 
functional may be constructed from this insight. We re- 
gard this as an initial step toward an accurate approxi- 
mation for /xc. Another obvious analytic tool would be 
the direct adiabatic decomposition of /xc in terms of A, 
which has proved so successful in understanding the hy- 
brids commonly used in ground-state calculations . 
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